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The invariant mass spectrum of the lepton pair ini-jnclusive semileptonic B —f Xu(.v decay 
yields a model independent determination of |K,j,|.u Unlike the lepton energy and hadronic 
invariant mass spectra, nonperturbative effects are only important in the resonance region, 
and play a parametrically suppressed role when AT /d(f is integrated over > (ms — mn)^, 
which is required to eliminate the B — > Xc£i> background. We discuss these backgrounds for 
slightly below (ms — mn)^, and point out that instead of > (m_B — mo)^ = 11-6 GeV^, 
the cut can be lowered to q^ > 10.5 GeV^. This is important experimentally, particularly 
when effects of a finite neutrino reconstruction resolution are included. 

A precise and model independent determination of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix element Vub is important for testing the Standard Model at B factories via the comparison 
of the angles and the sides of the unitarity triangle. At the present time the allowed range for 
sin 2/? in the SM is largely controlled by the model dependent theory error in \Vub\- 

If it were not for the huge background from decays to charm, it would be straightforward 
to determine \Vub\- Inclusive B decay rates can be computed model independently in a series 
in AqcD/n-b and as{mi,) using an operator product expansion (OPE)pEl and the result may 
schematically be written as 

dr=f^7^^^)x(i + A + Z%M + ... + ^!£(...) + 4(.. .) + ...). (1) 

V decay J [ rrib mf, vr vr^ J 

At leading order, the B meson decay rate is equal to the b quark decay rate. The leading 
nonperturbative corrections of order Aqq^/tti^ are characterized by two heavy quark effective 
theory (HQET) matrix elements, usually called Ai and A2. These matrix elements also occur in 
the expansion of the B and B* masses in powers of AqcD/^ife, 

T Ai + 3A2 7 Ax — A2 

rriB = rub + A h • • • , niB* = rub + A \- . . . . (2) 

Zrub Zrub 

Similar formulae hold for the D and D* masses. The parameters A and Ai are independent 
of the heavy b quark mass, while there is a weak logarithmic scale dependence in A2. The 
measured B* — B mass splitting fixes A2(mfe) = 0.12 GeV^, while A and Ai (or, equiva^ntly, a 
short distance b quark mass and Ai) may be determined from other physical quantities .B"EI Tims, 
a measurement of the total B — >■ Xuiu rate would provide a ~ 5% determination of 
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Figure 1: The shapes of the lepton energy and hadronic invariant mass spectra. The dashed curves are the 6 
quark decay results to 0{as), while the solid curves are obtained by smearing with the model distribution function 
f{k+) in Eq. (H). The unshaded side of the vertical lines indicate the region free from charm background. 



Unfortunately, the B X^H-i^ rate can only be measured imposing cuts on the phase space 
to eliminate the ~ 100 times larger B — s- Xciu background. The predictions of the OPE are 
only model independent for sufficiently inclusive observables, when hadronic final state with 

m\ » ^xAqcd > Aqcd (3) 

are allowed to contribute. Two kinematic regions for which the charm background is absent 
have received much attention: the large lepton enei:gY_region, Ef^ > {rn?^ — m'j^)/2mB, and the 
small hadronic invariant mass region, mx < mjjS^''^ However, in both of these regions of 
phase space the B — > Xuli? decay products are dominated by high energy, low invariant mass 
hadronic states, for which the inequality (^) is violated and the OPE breaks down. This occurs 
because the OPE includes the expansion parameter Ex-^qcb/^'x which becomes of order unity 
(mfeAqcD/w-c ^ 1 numerically) for Ex ~ rrii, and mx ~ rric- To predict the rates in these 
regions, the complete series in Ex^qcd/^x niust be resummed into a nonperturbative light- 
cone distribution function f{k^) for thefe nuark.El To leading order in l/rrife, the effects of the 
distribution function on various spectrat3'EJ may be included by replacing mi, by m^ = m^ + A;+ 
in the parton level spectrum, dPp, and integrating over the light-cone momentum 

dr = / dk+ f{k+) dPp , . (4) 

The situation is illustrated in Fig. ||, where we have plotted the lepton energy and hadronic 
invariant mass spectra in the parton model (dashed curveaLand smeared with a simple one- 
parameter model for the distribution function (solid curves 
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6(1 -x), ^ = X' A = 0.48GeV. (5) 



While it may be possible to extract f{k^) from the B Xgj photon spectrum,llj'llj unknown 
order Aqcd/"?-6 corrections are left over, limiting the accuracy with which \ Vub\ may be obtained. 

The situation is very different for the dilepton invariant mass spectrum. Decays with q'^ = 
{pi + Pp)'^ > {rriB — w-d)^ must arise from 6 — s- li transition. Such a cut forbids the hadronic 
final state from moving fast in the B rest frame, and simultaneously imposes mx < Tnjj and 
Ex < mo- Thus, the light-cone expansion which gives rise to the shape function is not relevant 
in this region of phase spaccOO This is also clear from Eq. (|^: the contribution of the Ai 
term to the decay rate, which is the first term in the shape function, is suppressed compared to 
the lowest order term in the OPE for any value of g^. The effect of smearing the q'^ spectrum 
with the model distribution function in Eq. (^) is illustrated in Fig. |2|. It is clearly a subleading 
effect. The improved behavior of the spectrum over the Ei and m'j^ spectra is also reflected in 
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Figure 2: The dilepton invariant mass spectrum. The notation is the same as in Fig. |l| 



the perturbation series. There are Sudakov double logarithms near the phase space boundaries 
in the Ei and spectra, whereas there are only single logarithms in the spectrum. 

The B — > Xu^u decay rate with lepton invariant mass above a given cutoff can be reliably 
computed working to a fixed order in the OPE (i.e., ignoring the light-cone distribution function), 



s2 _ „2/^2 /3_ _ 11 o^./Q r_ — n2 IT/.I2, 



where qr = q /ml, (3^ = 11 — 2nj/3, and Fq = Gp |T4bp m^/(192 vr^) is the tree level h ^ u 
decay rate. The ellipses in Eq. @ denote terms of order (Aogp/mf,)^ and order a1 terms not 
enhanced by Bn. The function X{q^) is known analytically ,t3 whereas Yia^) was computed 
numerically J13 The order l/m\ nonperturbative corrections are also known,c3 as are the leading 
logarithmic perturbative corrections proportional to a" log"'(mc/mb).E2l The matrix element of 
the kinetic energy operator, Ai, only enters the spectrum in a very simple form, because the 
unit operator and the kinetic energy operator are related by reparameterization invariance.H 

The relation between the total B — > Xuii? decay rate and \Vub\ is known at the 
levell'l^ 



'B{B XJiy)\g2^g2 i.6ps^ 
0.001 X F(g2) 



iKbl = (3.04 ±0.06 ±0.08) X lO'M — Ev^f^ — ^ > (7) 



where F{qQ) is the fraction of B ^ X^^lv events with > Qq, satisfying F{0) = 1. The 
errors explicitly shown in Eq. (0) are the estimates of the perturbative and nonperturbative 
uncertainties in the upsilon expansiona respectively. At the present time the biggest uncertainty 
is due to the error of a short distance b quark mass, whichever way it is definedJHll (This can be 
cast into an uncertainty in an appropriately defined A, or the nonperturbative contribution to the 
T(15') mass, etc.) By the time the spectrum in B -^JCu^^ is measured Jhis uncertainty should 



be reduced from extracting mb from the hadron massB or lepton energyS spectra in S — > Xcf-v, 
or from the photon energy spectrum! in i? ^ Xg^. The uncertainty in the perturbation theory 
calculation will be largely reduced by computing the full order correction in Eq. The 
largest "irreducible" uncertainty is from order Aq^j-j/m^ terms in the OPE, the estimated size 
of which is shown in Fig. ^, together with our central value for F(gg), as functions of q^. 

There is another advantage of the spectrum over the mx spectrum to measure \Vub\- In 
the variable mx, about 20% of the charm background is located right next to the 6 — s- u "signal 
region", mx < namely B — > Dii' at mx = mn- In the variable g^, the charm background 
just below q^ = {ms — mrif' comes from the lowest mass Xc states. Their q^ distributions 
are well understood based on heavy quark symmetry since this region corresponds to near 
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Figure 3: (a) The fraction of B ^ Xuiv events with > go, F(qQ), in the upsilon expansion. The dashed line 
indicates the lower cut = (m_B — mo)'^ — 11.6 GeV^, which corresponds to F = 0.178 ± 0.012. The shaded 
region is the estimated uncertainty due to Aqcd/ti^ terms; which is shown in (b) as a percentage of F{q^). 
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Figure 4: Charm backgrounds near q^ = (m_g — raoY ■ (Arbitrary units.) 



zero remil. Fig. |^ shows the B Dlu and B D*li> decay rates using the measured form 
factor^i (and \Vuh\ = 0.0035). The B XuC-v rate is the flat curve. Integrated over the region 
(f' > {iriB — mu*)'^ ~ 10.7GeV^, the uncertainty of the B ^ D background is smah due to its 
{w"^ — 1)^/^ suppression compared to the B — > X^iv signal. This uncertainty wih be further 
reduced in the near future. This increases the h ^ u region relevant for measuring \Vub\ by 
~ IGeV^. The B D* rate is only suppressed by {w'^ — 1)^/^ near zero recoil, and therefore 
it is more difficult to subtract it reliably from the b ^ u signal. The nonresonant Dtt final 
state contributes in the same region as .B — > D*, and it is reUably predicted to be small near 
maximal (zero recoil) based on chiral perturbation theory .E3 The D** states only contribute 
for q"^ < 9GeV^, and some aspects of their spectra are also known model independently .H 

Concerning experimental considerations, measuring the spectrum requires reconstruction 
of the neutrino four-momentum, just like measuring the hadronic invariant mass spectrum. A 
lepton energy cut may be required for this technique, however, the constraint q^ > {niB — mD)'^ 
automatically implies E£ > {tub — niD)'^ f^mB — 1-1 GeV in the B rest frame. Even if the Ei 
cut has to be slightly larger than this, the utility of our method will not be affected, but a 
calculation including the effects of arbitrary E£ and q^ cuts would be required. If experimental 
resolution on the reconstruction of the neutrino momentum necessitates a significantly larger 
cut than q^ = {uib — uibi)'^, then the uncertainties in the OPE calculation of F{qQ) increase. In 
this case, it may be possible to obtain useful model independent information on the q'^ spectrum 
in the region q"^ > 'm'^^2s) — 13.6 GeV^ from the q^ spectrum in the rare decay B —>■ Xs£~^i~, 
which may be measured in the upcoming Tevatron Run-II. 

In conclusion, we have shown that the q^ spectrum in inclusive semileptonic B — > X^lv decay 



gives a model independent determination of \Vub\ with small theoretical uncertainty. Nonpertur- 
bative effects are only important in the resonance region, and play a parametrically suppressed 
role when dT/dq'^ is integrated over q'^ > {tub — m^))^, which is required to eliminate the charm 
background. This is a qualitatively better situation than other extractions of \ Vub\ from inclusive 
charmless semileptonic B decay. 
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